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Abstract. For a fixed virtual scene (=collection of simplices) <S and given observer position 
p, liow many elements of S are weakly visible (i.e. not fully occluded by others) from p? The 
present work explores the trade-off between query time and preprocessing space for these 
quantities in 2D: exactly, in the approximate deterministic, and in the probabilistic sense. 
We deduce the existence of an 0{m^ /n^) space data structure for <S that, given p and time 
0{\ogn), allows to approximate the ratio of occluded segments up to arbitrary constant 
absolute error; here m denotes the size of the Visibility Graph — which may be quadratic, but 
typically is just linear in the size n of the scene S. On the other hand, we present a data 
structure constructible in 0(n ■ log(n) + ■ polylog(n)/^) preprocessing time and space 
with similar approximation properties and query time 0{£ ■ polylogn), where 1 < £ < n is 
an arbitrary parameter. We describe an implementation of this approach and demonstrate 
the practical benefit of the parameter £ to trade memory for query time in an empirical 
evaluation on three classes of benchmark scenes. 
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1 Motivation and Introduction 

Back in the early days of computer graphics, hidden surface removal (and visible surface calculation) 

was a serious computational problem: for a fixed virtual 3D scene and given observer position, 
(partition and) select those scene primitives which are (and are fully) visible to the observer. 
Because of its importance, this problem has received considerable scientific attention with many 
suggestions of deep both combinatorial and geometric algorithms for its efficient solution. The sit- 
uation changed entirely when the (rather unsophisticated) z-buffer algorithm became available in 
common consumer graphics cards: with direct hardware support and massive parallelism (one gate 
per pixel) , it easily outperforms software-based approaches with their (usually huge factors hidden 
in) asymptotic big-Oh running times |McKe87] . For a fixed resolution, the z-buffer can render 
scenes of N triangles on-line in time essentially linear in N with a small constant. However even 
this may be too slow in order to visualize virtual worlds consisting of several hundreds of millions 
of triangles at interactive frame rates. Computer graphics literature is filled with suggestions of 
how to circumvent this problem; for example by approximating (in some intuitive, informal sense) 
the observer's views. Here, the benefit of a new algorithm is traditionally demonstrated by eval- 
uating it, and comparing it to some previous 'standard' algorithm, on few 'standard' benchmark 
scenes and on selected hardware. We on the other hand are interested in algorithms with provable 
properties, and to this end restrict to 

1.1 Conservative Occlusion Culling 

Definition 1. Objects which are hidden to the observer behind (possibly a collection of) other 
objects may, but need not, be filtered from the stream sent to the rendering hardware, whereas any 
at least partially visible object must be visualized. 

Here, "conservative" refiects that the rendering algorithm must not affect the visual apprearance 
compared to the bruce-force approach of sending all objects to the hardware. Occlusion culling 
can speed up the visualization particularly of very large scenes (e.g. virtual worlds as in Second 
Life or World of Warcraft) where, composed from literally billions of triangles, typically 'just' some 
few millions are actually visible at any instant. Other scenes, or viewpoints within a scene, admit 
no sensible occlusion; for instance the leaves of a virtual forest naturally do not fully screen sight 
to the sun or to each other, similarly for CAD scenes of lattice or similar constructions. In such 
cases, spending computational efforts on occlusion culling is futile and actually bound to a net 
performance loss. Between those extremes, and particularly for an observer moving between oc- 
cluded and free parts of a large scene, the algorithmic overhead of more or less thoroughly filtering 
out hidden primitives generally trades off against the benefit in reduced rendering complexity. 
Put differently: Graphics hardware taking care of the visibility problem anyway opens the chance 
to hybridize with software performing either coarse (and quick) or careful (and slow) culling and 
leave the rest to the z-buffer. 

1.2 Adaptive Occlusion Culling 

It is our purpose is to explore this trade-off and to make algorithms adapt to each specific virtual 
scene and observer position in order to exploit it in a well-defined and predictable way. To this 
end we propose so-called visibility counts (the number of primitives weakly visible from a given 
observer position) as a quantitative measure of how densely occluded a rendering frame is, and 
whether and by how much occlusion culling therefore can be, expected to pay off. For technical 
reasons employed in Sections 12.61 and later, the formal notion slightly more generally captures the 
visibility of 'target' scenes through 'occluder' scenes: 

Definition 2 (visibility count). For a scene S — {Si, . . . ,5'„} of ^geometric primitives' Si C 
R'^, a subset of 'targets ' T C_ S, and an observer position p G M'', let 

ViS,p,T) :^ {T eT\3qeT -.yS e S\{T} : [p,q]° n S ^ID} 
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and denote by Y{S,p,T) := CardV(5,p, T) the number of objects in T weakly visible (i.e. not 
fully occluded) from p through S. Here, [p, q\° := {X ■ p + {1 — X) ■ q : < \ < 1} denotes the 
(relatively open) straight line segment connecting p and q. 

For scenes S and observer positions x with V{S,x,S) <C Card 5, occlusion culling is likely to pay 
off; whereas for Y{S,x,S) = Card 5 it is not. Quantitatively we have the following 

Hypothesis 3. Each culling algorithm A can be assigned a threshold function e [0, 1] such 

that, for scenes S and observer positions x with visibility ratios V(5, x,S)/ Card(5) (significantly) 
beyong 0j[{n), it yields a net rendering benefit and (significantly) below does not. 

1.3 Combinatorial Geometry and Randomized Computation 

Adaptivity constitutes an important issue in Computational Geometry; for instance in the context 
of Range Searching problems whose running time is preferably output sensitive, i.e. of the form 
©(/(n) + fc) where n denotes the overall number of objects and k those that are actually reported; 
compare |Chaz86] . 

Adaptivity is of course a big topic in computer graphics as well. However this entire field, 
driven by the impetus to quickly visualize (e.g. at 20fps) concrete scenes in newest interactive 
video games, generally focuses on innovative heuristics and techniques at a tremendous pace. We 
on the other hand are interested in algorithms with provable properties based on formal and sound 
analyses and in particular with respect to well-defined measures of adaptivity. This of course calls 
for an application of computational and combinatorial geometry [BKOS97|Edel87j . 

Paradigm 4 (Computational Geometry in Computer Graphics). For interactive visual- 
ization of very large virtual scenes of size N > 10^, algorithms must run in sublinear time 0{N"), 
a < 1, using preprocessed data structures of almost linear space 0(N^'^'^), e <C 1, provably! 

Here (time and) space complexity refers to the number of (operations on) unit-size real coordinates 
used (performed) by an algorithm — as opposed to, e.g., rationals of varying bitlengtl l* * *l . Also 
visibility is considered in the geometric sense (as opposed to e.g. pixel-based notions): point q 
is visible from observer position p if both can be connected by an ideal light ray (=straight-line 
segment not intersecting any other part of the scene), recall Definition [51 Our algorithm features 
a parameter l<£<n to trade preprocessing space for query time. 

Randomized algorithms are quite common in computer science for their efficiency and imple- 
mentation simplicity. They have also entered the field of computer graphics. Here these techniques 
are employed to render only a small random sample of the (typically very large) scene in such a 
way that it appears similar to the entire scene |WFP*01|KKF*04|WW*06] . Our goal, on the other 
hand, is to approximate the count of visible objects (Definition [2]) , not their appearance. 

1.4 Visibility 

Visibility comprises a highly active field of research, both heuristically and in the sound frame- 
work of computational geometry [CCSD03] . Particularly the latter has proven combinatorially 
and algorithmically non-trivial already in the plane [ORou87IGhos07] . Here the case of (simple) 
polygons is well studied |CAF07| : and so is point -point, point-segment, and segment -segment 
visibility for scenes 5 of n non-crossing line segments, captured e.g. in the Visibility Graph data 
structure [ GhMo91J . Its nodes correspond either to segments or to segment endpoints (or to both: 
a bipartite graph); and two nodes get joined by an edge if one can partly see the other. Weak 
segment-segment visibility for instance amounts to the 0{n'^) questions (namly for each pair of 
segments A and B) of whether there exist points a e A and b B such that a is visible from b. 

We, too, ask for weak segment visibility; however in our case the observer is not restricted to 
positions on segments of the scene but may move freely between them. For instance we shall want 
to efficiently calculate visibility counts for singleton targets V(5, x, {T}): 



* * * See however Item d) in Section [S] below 
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Problem 5. Fix a collection S of non-intersecting segments in the plane and one further segment 
T. Preprocess {S, T) into an almost linear (or merely worst-case subquadratic) size data structure 
such as to decide in sublinear time queries of the following type: 
Given x G , is T (partly) visible through S ? 

Sections 12. li 12. 2i and l2.3l recall two algorithms that meet either the space or the time requirement 
but not both. 

1.5 Overview 

An empirical verification of Hypothesis [3] in dimensions 2, 2^, and 3, is the subject of a separate 
work [FJZ09] . Our aim here is to explore the complexity of calculating the visibility counts, thus 
providing rendering algorithms with the information for deciding whether to cull or not. In view of 
the large virtual scenes and the high frame rates required by applications, we have to consider both 
computational resources, query time and preprocessing space, simultaneously. Section [2] focuses on 
the problem of calculating visibility counts exactly, mostly based on the Visibility Space Partition; 
our main result here is a preprocessing algorithm with output-sensitive running time. Section [3] 
weakenes the problem to approximate calculations: first showing the existence of a rather small 
data structure with logarithmic query time in Section [3.1l However this data structure seems hard 
to construct in reasonable time, therefore Sections I3.2f f consider approaches based on random 
sampling. Section [3] describes an implementation and evaluation of this algorithm. 

2 Exact Visibility Counting 

This section recalls combinatorial worst-case approaches for calculating visibility counts according 
to Definition [2l Many efficient algorithms are known for visibility reporting problems, that is for 
determining the view of an observer [Pocc90j ; however since reporting may involve output of linear 
size, such aproaches are generally inappropriate for our goal of counting in sw&linear time. On the 
other hand, logarithmic time becomes easily feasible when permitting quartic space in the worst- 
case based on the Visibility Space Partition (VSP). The main result of this section. Theorem [T3l 
yields an output-sensitive time algorithm for computing the VSP of a given set of line segments 
in the plane. 

2.1 Reverse Painter's Algorithm 

Prior to the hardware z- buffer, Painter's Algorithm was sometimes considered as a means to hidden 
surface elimination (at least in the 2D case): Draw all objects in back-to- front order, thus making 
closer ones paint over (and thus correctly cover) those further away. This of course relies on being 
able to efficiently find such an order: which is easily seen impossible in general unless we 'cut' 
some objects. Now two-dimensional BSP Trees provide a means to find such an order and a way to 
cut objects appropriately without increasing the overall size too much. We report from [BKOS97| 
Section 12]: 

Fact 6. Given a collection S of n non-crossing line segments in the plane, a BSP Tree of S can 
be constructed in time and space 0{n ■ logn). 

Now instead of drawing the cut segments in back-to- front order (relative to the observer) , feeding 
them into an Interval Tree in front-to-back order reveals exactly which of them are weakly visible 
and which not. Since insertion into an Interval Tree of size n takes time 0{n ■ logn) we conclude 
|TeSe91| : 

Lemma 7. Given a collection S of n non-crossing line segments in the plane and an observer 
position x G M^, Y{S,x,S) can be calculated in time 0{n ■ log^ n) and space 0{n ■ logn). 

Notice that preprocessing S into a BSP Tree accelerates the running time 'only' by a constant 
factor. 
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2.2 Rotational Sweep 

One can improve Lemma [7] by a logarithmic factor: 

Lemma 8. Given a collection S of n non-crossing line segments in the plane and an observer 
position X, \{S,x^S) can he calculated in time 0{n ■ logn) and space 0{n). 

Proof. Sketch First mark aU segments invisible. Then consider the 2n endpoints of 5 in angular 
order around x while keeping track of the order of the segments according to their proximity to 
the observer, the closest one thus being visible: whenever a new segments starts insert it into an 
appropriate data structure in time O(logn), whenever one ends remove it. Since the initial sorting 
also takes time 0{n ■ logn), we remain within the claimed bounds. □ 

Nevertheless the running time still fails to meet Paradigm U) Also, these approaches seem to offer 
no way to take advantage of a singleton target for the purpose of Problem [5l 

2.3 Visibility Space Partition 

Lemmas [7] and [8] work without any, and do not benefit asymptotically from, preprocessing of the 
fixed scene S. On the other hand by the so-called locus approach — storing all visibility counts in a 
Visibility Space Partition (VSP) — they can later be recovered in logarithmic running time [SchiOlj : 

Lemma 9. a) For a collection S of n non-crossing line segments in the plane, there exists a 
partition ofM.^ into 0{n'^) convex cells such that, for all observer positions x £ C within one 
cell C, V{S, -,5) is the same. 

b) The data structure indicated in a) and including for each cell its corresponding visibility count 
V(5, C,5) uses storage 0{n'^ ■ logn) and can he computed in time 0{n^ ■ logn). Then given 
an observer position x, its corresponding cell C , and the associated visibility count, can be 
identified in time ©(logn). 

c) When charging only real data and operations (more specifically: If an n-bit string is considered 
to occupy one memory cell and the union of two of them computable within one step), the 
above data structure including for each cell its visibility V{S,C,S) uses storage 0{n'^) and can 
be calculated in time 0{n^ ■ log^ n). 

d) Item a) extends to the case of [d — l)-simplices in d-dimensional space in that the number of 
convex cells with equivalent observer visibility can he hounded by 0{n)'^ . 

Proof, a) Draw lines through all (^2") pairs of the 2n segment endpoints. It is easy to see that, in 
order for a near segment to appear in sight, the observer has to cross one of these ©(n^) lines; 
compare Lemma [T51 below. Hence, within each of the 0{n'^) cells they induce, the subset of 
segments weakly visible remains the same; compare Figure [1] 

b) O(n^) lines induce an arrangement of overall complexity, and can be constructed in time 
©(n**) [BKOS971 Section 8.3]. The visibility number associated with each ceU is bounded by 
n and hence can be stored using ©(logn) bits; its calculation according Lemma [8] takes time 
0{n ■ logn) each. Finally, the planar subdivision induced by the 0{n*) edges of the arrange- 
ment can be turned into a data structure supporting point-location in ©(logn) ^BKOS97[ 
Theorem 6.8]. 

c) In the proof to b), constructing the arrangement (i.e. the planar partition into cells) and 
determining the visibility count of each cell were two separate steps which we now merge using 
divide- and-conquer: In the first phase calculate the VSP of the first two segments of S, then 
that of the next, and so on; in each VSP store, for each cell, the visibility vector, i.e. the 
0/1 bitstring recording which segments from S are visible (1) and which are not (0). In the 
next phase overlay the first two VSPs of two segments into one of the first four segments, and 
store for each refined cell the union of the 0/1 bitstrings: thus keeping track of its visibility; 
similarly for the next two VSPs of the next four segments. Then proceed to VSPs of eight 
segments each; and so on. We therefore have ©(logn) phases; and, according to [BKOS97| 
Theorem 2.6], the last (as well as each previous) phase takes time 0{n'^ ■ logn). 



6 M. Fischer, M. Hilbig, C. Jahn, F. Meyer auf der Heide, M. Ziegler 

d) Similarly to the proof of a), consider all dn vertices of the n simplices. Any d-tuple of them 
induces a hyperplane; and a change in sight requires the observer to cross some of these N :— 
i^f) ^ 0{nY hyperplanes. N hyperplanes in d-space induce an arrangement of complexity 
0{NY jEdel87) . □ 



Fig. 1. Visibility Space Partition of three segments: an observer in the dark gray area can see 
exactly one segment, in the light gray area exactly two, and otherwise all three of them. 



2.4 Size of Visibility Space Partitions 

Lemma in^+b) bounds the size of the VSP data structure by order n^. It turns out that this bound 
is sharp in the worst case — but not for many 'realistic' examples. This is due to many of the 0{n^) 
lines employed in the proof of Lemma [9^) inducing unnecessarily fine subdivisions of viewpoint 
space. In Figure [3^) for instance, the dotted parts are dispensible. 

In order to avoid trivialities, we want to restrict to nondegenerate segment configurations 
S. However this notion is subtle because the lines induced by S defining the VSP typically are 
degenerate: many (more than two) of them meet in one common (segment end) point. 

Definition 10. A family S of segments in the plane is nondegenerate if 

i) any two segments meet only in their common endpoints. 

ii) No three endpoints share a common line; 

Hi) Any two lines, defined by pairs of endpoints, do meet. 

We have already referred to (and implicitly employed in Lemma [2] a refinement of) the Visibility 
Space Partition; so here finally comes the formal 

Definition 11. For two non- degenerate collections S and T of segments in the plane, partition 
all viewpoints p G into classes having equal visibility V{S, p, T) . Moreover let VSP (5, T) denote 
the collection of connected components of these equivalence classes. The size of VSP is the number 
of line segments forming the boundaries of these components. 

Observe that VSP(5, T) indeed constitutes a planar subdivision: a coarsening of the 0{n'^) convex 
polygons induced by the arrangement of 0{n'^) lines from the proof of Lemma [9^). In fact a class 
of viewpoints of equal visibility can be disconnected and delimited by very many segments, hence 
merely counting the number of classes or cells does not refiect the combinatorial complexity. 
Lemma |9^) and Lemma [12] a) correspond to |Mato02[ Exercise 6.1.7]. 



Lemma 12. a) Even for a singleton target T , there exist a nondegenerate line segment configu- 
rations S such that VSP(5, {T}) has Qiji^) separate connected components. 
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b) To each n, there exists a nondegenerate configuration S of at least n segments admitting a 
convex planar subdivision of complexity 0{n) such that, from within each cell, the view to S 
is constant; i.e. VSP(5,iS) has linear size. 

c) The size of\SV[S,S) is at most quadratic in the size m of the Visibility Graph of S (recall 
Section \1.4-^ . 

d) A data structure as in Lemma\^ can be calculated in time 0{n ■ logn + • log^ n) and space 
0{m^). 

Since the Visibility Graph itself can have at most quadratically more edges than vertices, Item c) 
strengthens LemmalH^). Empirically we have found that a 'random' scene typically induces a VSP 
of roughly quadratic size. This agrees with a 'typical' scene to have a linear size Visibility Graph 
according to |ELPZ07| . 

Proof (Lemma\Wjj. 

a) Figure [2] is a small modification of |ORou87| FiG. 8.13]. The long bottom line segment T is 
visible from the upper half iff the observer can peep through two successive gaps simultaneously, 
i.e. from any position on the 0{n'^) stripes but not from the ellipses. When moving from an 
ellipse to another, T flashes into sight and is then hidden again. There are (9(n^) such ellipses. 
It is easy to see that this example is combinatorially stable under small perturbation and hence 
can be made non-degenerate. 

b) Consider Figured Within each cell C of the segment arrangement, all segments are always 
visible; hence C is also a cell of the VSP. And the exterior of S gives rise to another 9 VSP 
cells. 

c) Recall the proof of Lemma|9^); but throw in the observation that crossing the line g induced 
by two endpoints a and b (say, of segments 81,82 € S, respectively) does not induce a change 
in visibility if a and b are occluded from each other by some further segment 8 G S: compare 
Figure [31 Hence it suffices to consider at most as many lines as the the number m of edges 
in the Visibility Graph; and these induce an arrangement of at most quadratic complexity 

d) Determine, according to [GhMo9l| in time 0{n-\ogn + m) — or maybe more practically in time 
0{m ■ logn) |OvWe88] — the m lines mentioned in Item c); then proceed as in Lemma[3J;). □ 




2.5 Output-Sensitive VSP Calculation 

In view of the large variation of VSP sizes from order n to order according to Lemma 1121 the 
algorithms indicated in Lemma [OJd+c) for their calculation are reasonable only in case of large 
VSPs. We now present an output-sensitive improvement of Lemma [T2ji): 
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Fig. 3. Visibility does not change to an ob- 
server crossing the dotted (parts of) Unes Fig. 4. Example of a segment arrangement whose 
through pairs of segment endpoints. Visibility Space Partition has only linear size. 



Theorem 13. In 2D, the data structure of Lemma\9p,-c) can he computed in time 0{n'^ ■ logn + 
N ■ logn) in the sense of (the computational model referred to in) Lemma\3jp), where N denotes 
the combinatorial complexity ofYST'{S,S). 

Proof. We start as in the proof Lemma [9] with the order lines induced by all pairs of segment 
endpoints. Now the idea is to extend Lemma [T^t), namely to take into consideration only those 
parts the lines lines are cut into, which to cross actually changes the visibility. Indeed, these 
sub-lines constitute the boundaries of the cells of the VSP and therefore determine its complexity. 

i) Take such a line g passing through segment endpoints a and b of segments Si and 5*2 G S. To 
an observer crossing L_. g, the visibility can, but need not, change — and we want to determine 

if and where it does. First observe that the middle part (a, b) of g can be disposed off right 

away (unless a and b are endpoints of the same segment, but these 0{n) cases give rise to 
only 0{n'^) combinatorial complexity anyway) because crossing it never changes the visibility; 
compare Figured 

Now consider the two remaining unbounded rays of L- g, starting from a and starting 

from b. If, say, Ls intersects some other segment 5' € 5 in some point c, then traversing the 
part of beyond that point does not affect the visibility either, as ^2 is 'shielded' from sight 
by S anyway; again cf. Figure [31 So let (a, c) in this case, LL ■= Lg otherwise; and 

similarly for Lt. 

Now crossing at some point may or may not alter the visibility of (at least one of) 81,82, S 
(the latter being a segment 'opposite^ to b along g) but if it does so, then it does so at every 
point of Lg. Hence we will either keep the whole Lg, or drop it entirely; similarly for L'g. 

ii) Now since those two alternatives — namely keeping or dropping L^ — depend only on 81,82, 8, 
they can be distinguished in constant time. Moreover, after 0{n'^) preprocessing time and 
space for S, each L- g can be decomposed into the two parts L'~ and Lg as the result of a ray 
shooting query among S in logarithmic time; see e.g. jPoccQOl Theorem 3.2]. 

The line parts L' kept will in general intersect each other. So next cut them into non- 
intersecting maximal sub-segments. By the above observations, these constitute the boundaries 
of the VSP. And as a standard segment intersection problem, they can be determined in time 
0{n^ - log n + N); cf. e.g. [BKOS97, Section 2.5]. 

The resulting (sub-)segments give rise to a planar subdivision. For instance they cannot contain 
leaf (e.g. degree-1) vertices: circling around such a vertex one way would change the visibility 
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and the other way would not. Therefore the data structure admitting logarithmic-time point- 
location in the VSP can be calculated in space 0{N) and time 0{N ■ logn), recall [BKOS971 
Theorem 6.8]. 

iii) Determining the visibilities as in the proof of Lemma [lb) yields a factor n overhead; and 
the divide- and- conquer approach of Lemma [9t) seems inapplicable because of the correlations 
between segments in Step ii) , namely cutting off induced by Si , S2 at the first further 
segment S hit. On the other hand, each L'^ (and similarly for L'j.) by construction induces 
a definite change in visibility when crossed: we may presume this information to have been 
stored with L'^ at the beginning of Step ii). Hence we may start at one arbitrary cell of the 
arrangement, calculate its visibility according to Lemma [71 and then traverse the rest of the 
arrangement cell by cell while keeping track of the visibility changes induced by (and stored 
with) each cell boundary. □ 



2.6 Visibility of One Single Target: Trading Time for Space 

The query time obtained in Lemma [9] is very fast: logarithmic (i.e. optimally) where, according 
to Paradigm m sublinear suffices. Quite intuitively it should be possible to reduce the memory 
consumption at the expense of increasing the time bound. We achieve this for the case of one 
target, that is the decision version of visibility x > V(5, x, {T}) € {0, 1}: 

Theorem 14. For each 1 < ( < n, Problem\B[ can be solved, after 0{n^ ■ log^ n/i) time and space 
0{n'^/£) preprocessing, within query time 0{£ ■ logrt). 

Such a trade-off result from time to space has become famous in the general context of structural 
complexity [HPV77j . Note that, obeying sublinear time, we can get arbitrarily close to cubic 
space — yet remain far from the joint resources consumption of an interval tree (Section l2.1|) . But 
first comes the already announced 

Lemma 15. Fix a collection 51+) {T} ofn+1 non-crossing segments in the plane. Let Li, . . . , 
denote the k = (^"'2'^^ lines induced by the pairs of endpoints of segments in S U {T}. For an 
observer moving in the plane, the weak visibility of T can change only as she crosses 

— either one of the lines Li intersecting T 

— or someline supporting a segment S G S; 

compare Figure\^ 

Proof. Standard continuity argument: Let p denote the observer's position and suppose point 
X £ T is visible, i.e. the segment [p, x\ does not intersect S € S. Now move p until x is just about 
to become hidden behind S" £ 5. Then start moving x onT such as to remain visible. Keep moving 
p and adjusting x: this is possible (at least) as long as the line through p and x avoids all endpoints 
of5U{r}. □ 



Proof ( Theorem \ 14^ . Consider, as in the proof of Lemma [9l the 0{n'^) lines induced by pairs of 
segment endpoints of S. Consider the intersections of these lines with T (if any). Partition T into 
0{£) sub-segments Ti, . . . ,Ti, each intersecting 0{n'^/£) of the above lines. For each piece T^, take 
the arrangement Ai of size 0((n^/£-|-n)^) induced by those lines intersecting Ti, and all 0{n) lines 
through one endpoint of Ti and one of some S (Iz S, and all 0(n) lines supporting segments from S. 
By Lemma [T51 within each cell C of Ai, the weak visibility of Ti is constant (either yes or no) and 
can be stored with C: Doing so for each Ai {I < i < £ < n) and each of the Oin^ !£?' -I- r? I£ -I- r?') 
cells C of A^i uses memory of order 0{n'^ /£ + + n'^£) — 0(n*/£) as claimed; and corresponding 
time according to Lemma [9};). 

Then, given a query point p e R'^, locating^ in each arrangement Ai takes total time O(^-logn); 
and yields the answer to whether Ti is weakly visible from p or not. Now T itself is of course visible 
iff some Ti is: a disjunction computable in another 0{£) steps. □ 



10 M. Fischer, M. Hilbig, C. Jahn, F. Meyer auf der Heide, M. Ziegler 




We even can combine Lemma [T2H ) with Theorem [14] to obtain 

Scholium 16. For each 1 < d- < n, Problem\^ can be solved, after preprocessing S in 0{n- log n + 
m? ■ \o^n/i) time into an 0{w? jt) size data structure, within query time 0{£ ■ logn) where m 
denotes the size (number of edges) of the Visibility Graph ofS. 

Proof. Instead of considering, and partitioning into i groups, all 0{n^) lines induced by pairs of 
segment endpoints, do so only for the 0{m) lines induced by pairs segment cndpoints visible to 
each other. □ 



3 Approximate Visibility Counting 

Lacking deterministic exact algorithms for calculating visibility counts satisfying both time and 
space requirements, we now resort to approximations: of Y{S,x,S) up to prescribable absolute 
error A: G N or, equivalently, of the visibility ratio V(5, 5)/ Card(5) up to absolute error e — 
k/ Card(5); recall Hypothesis [3l 

Remark 17. Relative errors make no sense as there is always a viewpoint x with V(iS, x,S) — 1. 

Corollarv[27l the main result of this section, presents a randomized approximation within sublinear 
time using almost cubic space in the worst-case and almost linear space in the 'typical' one. 



3.1 Deterministic Approach: Relaxed VSPs 

Visibility space partitions, and the algorithms based upon them, are so memory expensive because 
they discriminate (i.e. introduces separate arrangement cells for) observer positions whose visibility 
differs by as little as one; recall Definition [TlJ It seems that considerably more (time and) space 
efficient algorithms may be feasible by partitioning observer space into (or merely covering it by) 
more coarse classes: 
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Definition 18. Fix /c €E N and collections S and T of non-intersecting segments in the plane. 
Some covering {Ci, . . . , C/} ofM? is called a k-relaxed VSP of{S,T) if 

\/l<i<I\/p,qeCr- yiS,p, T) - V(5, q,T) < k . 

In the sequel we shall restrict to k-relaxed VSPs which constitute planar subdivisions (i.e. each Ci 
being a simple polygon); and refer to their size in the sense of Definition \ll[ 

Indeed, such VSPs allow for locating a given observer position x in logarithmic time to yield 
a cell Ci 3 X which, during preprocessing, had been assigned a value V(5, q, T) approximating 
V(5,x, T) up to absolute error at most k. 

Example 19. For CardT < /c, the trivial planar subdivision {M?} is a fc-relaxed VSP of {S,T). 
In particular the quartic lower size bound of Lemma 112b ) applies only to 0-relaxed VSPs but 
breaks down for fc > 1. 

This example suggests that much smaller (e.g. worst-case quadratic) sizes might become feasible 
when considering fc-relaxed VSPs for, say, fc « ^yn or even k ~ n/logn. Indeed we have the 
following lower and upper bounds: 

Proposition 20. a) For each n, k there exists a non-degenerate family S of segments in the 
plane such that any k-relaxed VSP has size at least n ■ \_(n — l)/(fc -I- 1)J . 

b) There also exist such families such that any k-relaxed VSP has size at least f2(n*/fc^). 

c) Let S be a non- degenerate family of n segments in the plane and N the size of its VSP. Then 
there exists a k-relaxed VSP of size \_N/(k -\- 1)J . 

d) There also exists a k-relaxed VSP of size 0{m'^ /k'^), where m denotes the size of the Visibility 
Graph of S. 

e) In fixed dimension d, ( Definition \18\ and) Item c) generalizes to k-relaxed VSPs of size 

Recall that N < m? < n^, thus leaving quadratic gap between a) and d) for k small; and between 
b) and d) for k large. Item c) succeeds over d) in cases where N asymptotically does not exceed 
m^/fc. 

Proof, a) Consider Figure [6] with n = 13 segments which obviously generalizes to arbitrary n. 
Moving from segment a along the arrow, each time crossing a dashed line amounts to an 
increase in visibility from {a} via {a,b}, {a, 6, c} and so on up to entire S. Hence to obtain 
cells of viewpoints with visibility varying by at most fc, we must keep at least every (fc -I- l)-st 
dashed border, that is \_{n — l)/{k-\- 1)J out of n — 1. By symmetry, the same argument applies 
when moving from segment b along the arrow, or from any other segment. 

b) A closer look at Figure [6] reveals it to induce a VSP of size /2(n^): For any segment, there 
are linearly many separate cones from which it can be seen through the gaps between the 
other segments; hence we have quadratically many cones, of which almost any two intersect; 
compare Figure [71 particularly its right part. 

Now in order to argue about fc-relaxed VSPs, replace in Figure [6] each single segment by fc 
scaled and shifted copies as indicated in the left of Figure [Tj that is we now have a scene of size 
N = n-k. Observe that, when entering and passing through a cone of visibility, the number of 
segments visible increases from its original value to an additional fc (drawn in levels of gray) . 
The visibility number thus varies by fc -I- 1, requiring a fc-relaxed VSP to subdivide the cone; 
indeed the entire cone! By the above considerations, these necessary boundaries induce an 
arrangement of complexity fl{n^) which, expressed in the size N of the scene, is n{N'^/k'^) as 
claimed. 

c) Classify the cells of VSP(5, S) according to their visibility count; and let Ni denote the number 
of boundary segments of VSP (5, 5) separating a cell with visibility count i from one with 
visibility count i + 1. Since any boundary segment does so for some i = l,...,n— 1, 

N = Ni+N2 + --- + Nn-1 

= (iVi + Nk+2 + A^2(fe+l) + l + ...) + (A^2 + Nk+3 + iV2(fe+i)+2 + . . .) 

+ {Ns + Nk+4 + A^2(fc+l)+4 + ...) + ••• + {Nk+1 + N2(k+1) + A^3(fc+1) + ■ • •) • 
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Fig. 6. Crossing each dashed hne increases visibiUty by one segment; hence keeping only every 
fc-th leads to visibility count variations within a cell of at least k. 

Hence by pigeonhole principle there exists 1 < k < fc + 1 such that [N/{k + 1)J > + 
Nk+i+K + A^2(fc+i)+K + • • • Now keep from VSP(5,5) exactly those boundary segments that 
either separate cells with visibility count k from ones with visibility count k + 1, or cells with 
visibility count k + l + k from ones with visibility count k + 2 + n, or cells with visibility count 
2(fc + 1) + K from visibility count 2(fc + 1) + k + 1 and so on. 

By the above considerations, this planar subdivision has complexity at most [N/ {k + 1)J . And 
by construction, it joins cells having visibility counts j ■ (k + 1) + k + 1, j • (k + 1) + k + 2, 
j ■ {k + 1) + K + 3, . . . , j ■ (k + 1) + {k + 1) + k; but maintains their separation from cells with 
visibility count (j + 1) • (fc+ 1) + k+ 1, because those boundary segments are precisely the ones 
deliberately kept. Hence the joined coarsened (i.e. super-) cells indeed contain only viewpoints 
with visibility count differing by at most k. 

d) Recall the collection £ of m lines used in the proof of LemmafT^). Now let r := m/ (fc— 1), k > 
1, and apply the Cutting Lemma of Chazelle, Friedman, and Matousek [Mato021 Lemma 4.5.3]: 

There exists a subdivision of the plane into C(r^) = 0{m'^/k^) generalized (i.e. not 
necessarily closed) triangles Ai such that the interior of each Ai is intersected by at 
most m/r = fc — 1 lines from £. 

Since visibility changes occur only when crossing lines in C, the latter means that the visibility 

counts within each Ai differ by at most fc. 

e) Instead of 0{r^) generalized triangles as in d), now employ simplicial cuttings of size 0{r'^) 
according to [Mato02[ Theorem 6.5.3]. □ 

We remark that Item c) considers a certain sub-arrangement of the 0-relaxed VSP, whereas the 
planar subdivision due to Item d) uses cell boundaries not necessarily belonging to the VSP. Also, 
the size 0{r'^) of the cuttings employed in Items d) and e) is known to be optimal in general; 
but this optimality does not necessarily carry over to our application in visibility, recall the gaps 
between lower and upper bounds in Items a) to d). 
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Fig. 7. Illustrating the proof of Proposition [20b) . 



Substituting k ^ S ■ n yields 

Corollary 21. For each collection S of non-degenerate segments in the plane and fc G N, there 
exists a data structure of size 0[min{N/{S ■ n),m'^/{S ■ n)^}) that allows to approximate, given 
X £ R^, the visibility ratio Y{S,x,S)l Gaid{S) up to absolute error S > in time 0{\ogn). 

Notice that we only claim the existence of such small data structures. In order to construct them, 
the proofs of Proposition [20b ) and d) both proceed by first calculating the 0-relaxed VSP and 
then coarsening it. Specifically for Proposition [20tl) . a Triangular Cutting can be obtained in time 
0{m ■ r ■ poly log) AgarQOJ ; but determining the visibility count for each triangle Ai costs 0{n ■ 
log n) according to Lemma HI or can be taken from the 0-relaxed VSP. The preprocessing time 
for Proposition [^Ob+d') thus is, up to polylogarithmic factors, that of Theorem [T31 i.e. roughly 
0{N): independent of, and not taking advantage of large values of, k. For the configuration from 
Lemma [T2k) for instance (recall Figure|2]), this results in a preprocessing time of order although 
the resulting 1-relaxed VSP has only size 0{n^). Alternatively, apply Lemma|8]to (one point from) 
each triangle Ai to obtain a running time of roughly output xn, that is still off optimal by one 
order of magnitude. And finally, the asymptotically 'small' size and time for calculating triangular 
cuttings hide in the big-Oh notation some large constants which are believed to prevent practical 
applicability. 



3.2 Random Sampling 

Both size and query time of the data structure due to CoroUarv 1211 are rather low; but because of 
the infavourable preprocessing time and hidden big-Oh overhead indicated above, we now proceed 
to random sampling, based on a rather simple generic algorithm: 

Algorithm 22. 

i) Guess a sample target T C S of size m. 

ii) Calculate the count Y{S,x,T) of objects in T visible through S. 
Hi) Return the ratio Y{S,x,T)/ Card(T); 

iv) and hope that it does not deviate too much from the 'true' value Y{S,x,S)/ Card(5). 
Item iv) is justified by the following 

Lemma 23. Fix x G M."^ and S > 0, then choose T C S as m independent identically distributed 
random draws from S. It holds 



Probr 



V{S,x,T)/m - V{S,x,S)/n\ > 6 



< 2-e'^'^-^^ 
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In other words: In Algorithm [21] taking m (quadratic in the aimed absolute accuracy d but) 
constant with respect to the scene size n suffices to achieve the desired approximation with constant 
probabiUty; shghtly increasing it further amplifies exponentially the chance for success. 

Remark 24- a) It is easy to see that a fixed relative accuracy can be attained, for Y{S, x, S)/n ^ 

0, only by samples of size m — > n: If only one segment is visible, it must get sampled to be 
detected. 

b) Also the visibility of the sample is crucially to be considered with respect to the entire scene, 

1. e. V(5,f,T) rather than V(T,f,T). 

Proof (Lemma \ 2S]} . It is well-known |MoRa95| AlSp00' that a sum X :— X^IHi "'^i '^^ indepen- 
dent {0,1} trials Xi,...,Xm satisfies the Chernoff-HoefFding Bound Prob [jAT/m — /i| > 6] < 
2 ■ exp(— 2m • 5^) where /i denotes the expectation of Xi. In our case, let Xi denote the event that 
the i-th draw St € S is visible from x through S. This happens with probability fi = V(5, x, S)/n, 
hence X = Y{S,x,T). □ 



3.3 The VC-Dimension of Visibility 

Note that the random experiment T and the probability analysis of its properties in Lemma [23] 
holds for each x but not uniformly in x. This means for our purpose to re-sample T C 5 at 
every frame. On the other hand, the above considerations have not exploited any geometry. An 
important connection between combinatorial sampling and geometric properties is captured by 
the Vapnik-Chervonenkis Dimension |AlSpOO| : 

Fact 25. Let X be a set and TZ a collection of subsets R C X. Denote by 

d := VCdim(X,7^) max { Card y |y C AT, {y n i? : i? G 7^} = 2^} (1) 
the VC-Dimension of{X,TZ). 

a) ForY CX,n^ Card(y), Card{y nR: ReU} < J^Lo id) ^ 

h) Let Y C X be random of Card(y) > max {A/S ■ log |, 8d/5 • log ^} . Then with probability at 
least 1 - p, it holds for each R G TZ: Card(A n R) > S ■ Card(A) ^ YDR^^d. 

c) Let Y <Z X be random o/ Card(y) > Qi^{d ■ log ^ -I- logi)/(5^). Then with probability at least 
l^pit holds for each RgTZ: \ Card(A nR)/ Card(A) - Card(y n i?)/ Card(y) \ < 6. 

Lemma 26. Fix a collection S of n non-crossing {d — l)-simplices in M.^ . 

a) Define X := S and TZ := {V(5,f,5) : x G R^}. Then VCdim(A,7^) < . (logn -f 0(1)). 

b) A random subset T CI S of cardinality m > i7(d^ • log(n) • log{d ■ \ogn/ 5) / 5) satisfies with 
constant (and easily amplifiable) probability that, whenever at least a 5- fraction of the simplices 
of S are visible from x G M."^ through S , then so is some simplex of T . 

c) A random subsetT C S of cardinality m > i7(d^ •log(n) •log((i-logn/(5)/(5^) satisfies with con- 
stant (yet easily amplifiable) probability thatY{S,x,T)/ Gaid{T) deviates from Y{S,x,S)/ Card(5) 
absolutely by no more than 6. 

d) The bound obtained in a) is asymptotically optimal with respect to n: In there exist non- 
degenerate collections S — X of n line segments such thatYGAim{X,TZ) > logn — ©(loglogn). 

Proof, a) Lemma [9H) imphes Card7^ < 0{n)'^\ Hence, in Equation 2^ ^ {Y n R : R e TZ} 

requires 2^'''"^^ < Card7^ and therefore Cardy < • (logn + 0(1)). 
b) and c) follow by plugging Item a) into Fact I^Sb-fc) . 

d) Figure [8] below obviously extends to the construction of k segments of which, using 2*^ • k/2 
additional segments as 'shields', each subset appears as a visible set V{S,x,S) for some x. 
That is a scene S of size n = k-\-2^ ■k/2 containing a subset Y of size k = logn — C'(loglogn) 
as in Equation (|T|). □ 
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c^:> ' tt^ 'Jft^ ' 

Fig. 8. 3 segments of which, shielded by 12 further segments, each of its 8 subsets appears as 
visible set. 



3.4 Main Result 

Lemma EBb ) enhances Lemma [521 The latter is concerned with the probabiUty of a constant-size 
sample T to be representative (i.e. to approximate the visibility ratio) with respect to a fixed 
viewpoint — i.e. our application would (have to) re-sample in each frame! The former lemma on 
the other hand asserts that a polylogarithmic-size sample, drawn once and for all, be suitable with 
respect to all viewpoints! In particular we may preprocess the visibihty of each T € T separately 
according to Theorem [14] and obtain, employing Scholium 1161 

Corollary 27. Given < S < 1, a collection S of n non-crossing segments in the plane (d — 2), 
and 1 < £ < n < m where m denotes the size of the Visibility Graph of S. Then a randomized 
algorithm can preprocess S within time 0{n ■ \ogn + nn? ■ polylogn • log j-/(i? ■ (5^)) and space 
0{w? ■ polylogn • \ogj/{(. -5^)) into a data structure having with high probability the following 
property: Given x G M.'^ , one can approximate the visibility ratio Y{S, x, S) / Card(iS) up to absolute 
error at most 5 in time 0{£ ■ polylogn • log \/ 5'^') ■ 

Again, note the tradc-off between space and query time gauged by the parameter £. And, remem- 
bering the paragraph following Lemma 1 121 m is 'typically' linear in n; hence choosing (. — n^~'^, 
the space can be made arbitrarily close to linear while maintaining sublinear query time, thus 
complying with Paradigm [¥] 



4 Empirical Evaluation 

The present section demonstrates the practical applicability of the algorithm underlying Corol- 
lary [57) It is (not trivial but neither) too hard to implement, constants hidden in big-Oh notation 
are modest, and query time can indeed be traded for memory. 

Measurements were obtained on an Intel® CPU 6700 running at 2.66GIIz under 

openSUSE 11.0 equipped with 4GB of RAM. The implementation is written in Java version 6 
update 11. Calculations on coordinates use exact rational arithmetic based on BigIntegers. 

4.1 Benchmark Scenes 

We consider three kinds of 'virtual scenes' in 2D, that is collections of non-intersecting line seg- 
ments, compare Figure [9] 

A) Sparse scenes representing forest-like virtual environments with long-range visibility; 
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1 1 



Fig. 9. Scene types A, B, and C. 



B) cellular scenes representing architectual virtual environments with visibility essentially limited 
to the room the observer is presently in; 

C) and an intermediate of both. 

As indicated by the above classification, these scenes contain some regularity. More precisely, their 
respective visibility ratios obey qualitative deterministic laws, see Figure ITUk ) . On the other hand 
these scenes are constructed using some random process, which means instances can be made up 
of any desired size n. 

Specifically all scenes arise from throwing into each square of an x ^yn grid one randomly 
oriented segment. For Scene A, these segments are then shrunk by a factor = ^/n to yield an average 
visibility count proportional to n, see Figure ITOk ). For Scene B, each segment sequentially is grown 
as to just touch some other one: remember we want to comply with Definition II 01 as expected (and 
corresponding to Lemma [T2b ) this results in constant visibility counts. Scene C finally arises from 
Scene B by shrinking the segments again; here the visibility count grows roughly proportional to 
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Fig. 10. a) Maximum/average/minimum visibility counts for Scenes A to C. 
b) Variance of the output of the randomized algorithm. 



Figure ITOb) indicates the quality of approximation attained by CoroUarv [27l Recall that the 
preprocessing step randomly selects k elements of the scene as targets; and the proof shows that 
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asymptotically this sample is 'representative' for the entire scene with high probability. Our im- 
plementation chose fc = 10 • log^(n) which turned out to yield a practically good approximation 
indeed. More precisely, Figure [TOb) displays the la confidence interval estimates for scenes A) and 
C) from various viewpoints, normalized to (i.e. after subtracting the) mean 0. 



4.2 Memory Consumption versus Query Time 

Preprocessing space and query time are two major resource contraints for many applications such 
as the one we aim at. We have thus performed extensive measurements of these quantities for 
the scene types A) to C) mentioned above. It turns out that for A) our data structure takes 
roughly linear space and for B) roughly quadratic one, whereas for C) it grows strictly stronger; 
see Figure ITTk) . Here we refer to the setting £ = 1. For scene C) we have additionally employed 
the trade-off featured by Corollary [27] to reduce the memory consumption at the expense of query 



time; specifically, scene E) means scene C) with £ and scene D) refers to £ 

turns out that the latter effectively reduces the size to quadratic, see Figure ITTk). 
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Fig. 11. Log-log plot of a) preprocessing memory and of b) query time for scenes A) to E). 



On the other hand, the saved memory is paid for by an increase in query time; cf. Figure [TTb l. 
Indeed at scenes of size n ^ 20, 000 each query is estimated to take about 50ms, that is as long 
as one frame may last at an interactive rate of 20fps. Whereas scene E), that is scene C) with 



i^^'^ instead off 



.1/2 



, is estimated to still remain far below this limit for much, much larger 



scenes n; not to mention scenes A) to C), i.e. with £ = 1. 



4.3 Conclusion 

Our benchmarks range up to n « 8, 000 when the data structure hit an overall memory limit of 
16GB. This may first seem to fall far short of the original sizes aimed at in Section [TTT] On the 
other hand, 

— the measurements obtained turn out to depend smoothly on n and thus give a sufficient 
indication of, and permit to extrapolate with convincing significance, the behavior on larger 
scenes. 

— By proceeding from single geometric simplicies to entire 'virtual objects' (like e.g. a house or 
a car) as rendering primitives, one can in practice easily save a factor of 100 or 1000. 

— Temporal and spacial coherence of an observer moving within a virtual scene suggests that 
visibility counting queries need not be performed in each frame separately. Moreover our CPU- 
based algorithm can be run concurrently to the graphics processing unit (GPU). These two 
improvements in running time can then be traded for an additional saving in memory. 
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— In order to attain and access the 16GB mentioned above, we employed secondary storage (a 
harddisk): with an unfavorable increase in preprocessing time but suprisingly little effects on 
the query time. 

These observations suggest that our algorithm's practicality can be extended to n much larger 
than the above 8, 000; yet doing so is beyond the purpose of the present work. 

It is thus fair to claim as main benefit of our contribution in Corollary 1271 an (as opposed to 
Corollary 1211) practically relevant approach to approximate visibility counting based on the abil- 
ity to trade (otherwise prohibitive quartic) preprocessing space for (otherwise almost neglectible 
logarithmic) query time. 

5 Perspectives 

a) We have treated the observer's position x as an input newly given from scratch for each frame. 
In practice however p is more likely to move continuously and with bounded velocity through 
the scene. This should be exploited algorithmically, e.g. in form of a visibility count maintenance 
problem [Pocc90l Section 3.2]. 

b) How does Theorem [Ml extend from 2D to 3D, what is the typical size of a 3D VSP? 

c) The quartic worst-case size of 2D VSPs (and quadratic typical yet even of order for 3D) 
arises from visibility considered with respect to perspective projections; whereas for ortho- 
graphic projections, it drops to 0{n^) (in 2D; in 3D: order n^) [SchiOlj . 

d) The counterexamples in Lemma [T2b) and Lemma [26H) and also Proposition[20k) employ (after 
scaling the entire scene to unit size) very short and/or very close segments. We wonder if such 
worst cases can be avoided in the bit cost model, i.e. with respect to n denoting the total 
binary length of the scene description on an integer grid. 

5.1 Remarks on Lower Bounds 

We have presented various data structures and algorithms for visibility counting, trading prepro- 
cessing space for query time. It would be most interesting to complement these results by cor- 
responding lower bounds of the form: preprocessing space s requires, in some appropriate model 
of computation, query time at least i7(/(s)). Unfortunately techniques for Friedman's Arithmetic 
Model, which have proven so very successful for range query problems [ChazQOj, do not apply 
to the non- faithful semigroup weights of geometric counting problems, not to mention decision 
Problem[51 and approximate, rather than exact, counting makes proofs even more complicated. 

On the other hand, we do have some lower bounds: namely on the sizes of VSPs and relaxed 
VSPs in Lemma [12k) and Proposition [20k-l-b') . These immediately translate to lower bounds on 
the sizes of Linear Branching Programs, based on the observation that any such program needs a 
different leaf for each different convex cell of inputs leading to the same output value, compare 
|DFZ02| . 

But then again this seemingly natural model of computation is put into question when con- 
sidering the algorithms from Sections 12.11 and 12.21 Both have, in spite of Lemma [T2k). merely 
(weakly) linear size. Superficially, Lemma [8] seems to employ transcendental functions for angle 
calculations; but these can be avoided by comparing slopes instead of angles — however employing 
divisions. These, again, can be replaced: yet multiplications do remain and make the algorithm 
inherently nonlinear a branching program. 

5.2 Visibility in Dimensions > 2 

We had deliberately restricted to the planar case of line segments. Many virtual scenes in inter- 
active walkthrough applications can be described as 2i-dimensional: buildings of various heights 
yet rooted on a common plane. 

But how about the full 3D case? Here we observe a quadratic 'almost' lower bound on the 
joint running time of preprocessing and querying for the 3D counterpart to Problem [5] To this 
end recall the following famous 
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Problem 28 (3SUM). Given an n-element subset S or M, do there exist a, 6, c £ S such that 

a + 6 + c = ? 

It admits an easy algebraic C'(n^)-time algorithm but is not known solvable in subquadratic time. 
Similar to Boolean Satisfiability (SAT) and the theory of A/'T'-completeness, 3SUM has led to a 
rich family of problems mutually reducible one to another in softly linear time 0{n ■ poly log n) 
and hence called ?>S\]M- complete; for example it holds [GaOv95. Section 6.1]: 

Fact 29. Given a collection S of opaque horizontal triangles in space, one further horizontal 
triangle T, and a viewpoint p e M"^. The question of whether some point of T is visible from p 
through S (called Visible-Triangle) is SSUM-complete. 

In particular there is no 3D counterpart to the Interval Tree solving the corresponding 2D problem 
in time 0{n ■ logn), recall Section [2T] 
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